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C Additional Focus (Porous Wick) 


While forced convection is a very effective method of cooling MCM's in terrestrial 
applications, it has a few drawbacks for space-based systems. The pump required to 
force the fluid through the heat exchanger adds weight to the system, and must also be 
electrically powered. 

Loop heat pipes are passive devices, which do not require a pump. Also, since they 
take advantage of the latent energy of the working fluid, they are able to dissipate a very 
large heat flux. Unlike capillary pumped loops (CPL's) they do not have evaporator start- 
up problems. The following description of how a loop heat pipe can be utilized to cool 
high power density MCM's is taken from a proposal [4] written by the Center for Space 
Power at Texas A&M University (Best, 1997). 
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Fig. 4: Loop heat pipe schematic 


“Figure (4) schematically shows a loop heat pipe composed of a compensation 
chamber an evaporative structure, capillary transport loops, and a capillary condenser 
Thermal energy from an MCM is transported to the evaporator wick by conduction and 
radiation. Evaporation from the wick produces vapor which travels through the vapor 
transport line to the condenser. The condenser transfers energy out of the working fluid 
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to a lower temperature external sink, converting the vapor to single phase liquid. The 
liquid returns to the evaporator through the liquid transport line which passes through the 
compensation chamber. The compensation chamber acts as an inventory controller for 
the system.” 

Currently one of the major limitations of a loop heat pipes ability to dissipate large 
heat fluxes is the porous wick. Conventional sintered particle wicks have a distribution 
of pore size and spacing which limit wick performance and hence maximum heat load for 
the loop heat pipe. Since the radii of curvature of the miniscus depends on pore size, the 
least favorable pore size determines the performance of the wick. 

Several years ago it was found that under the proper conditions an HF-based 
electrolytic fluid will etch a porous "sponge-like" structure in silicon, where the walls, 
re main single crystalline silicon, as shown by x-ray analysis. The porosity will show an 
anisotropic directionality lying along the predominant current flow direction, 
perpendicular to the semiconductor surface. In recent years, several research groups have 
pursued the electroluminescense properties where the surface will continue to glow after 
being initially exposed to intense light. However, the controllability of the pore sizes, the 
porosity (volume occupied by etched-out open space) and the thermal properties are of 

interest in this work. 


Also etching activity can be lithographically patterned as with integrated circuits, 
and supplemental control can be had by adjusting the molarity of the etchant, the driving 
current, the local crystal doping (more activity in p+ regions) and even the intensity ot 
directed light (as with a laser). Also, one can grow silicon dioxide on the surfaces ot the 
pore walls under a high temperature oxide atmosphere (as we do in microelectronic 
frabrication), thus since oxide adds to the wall thickness we can fine-time the pore size^ 
Of course the oxide thermal properties are different (lower) than that of silicon, thus w 
can adjust the thermal properties after fabrication. 

Micro-fabricating a porous wick in silicon allows one to create pores of uniform 
size and spacing. Also, it is possible to vary the diameter of individual pores across 
thickness of the wick. This may greatly increase the maximum heat load for a layered 
heat pipe by preventing the miniscus from receding all of the way into the liquid. The 
goal of this part of our research is to create pores with diameter of order 1 micron an 
thereby greatly improve the thermal performance of loop heat pipes. 
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in. PAST YEAR IN REV JEW (Accomplishments) 


B. Porous Wick Structure 


1. Microfabrication of Porous Silicon: Background Information and Technique 

In order to improve the performance of loop heat pipes it is necessary to fabricate a 
wick with a pore size of approximately 1 micron. Traditionally, porous silicon etching 
results in a sponge-like structure with pore sizes distributed from nanometers to hundreds 
of nanometers. 

On the other hand, macro porous silicon can create arrays of straight and uniform 
holes through a flat silicon wafer. This type of porous structure is expected to yield better 
performance than the usual sponge-like ones. The uniform cylindrical shape also makes 
modeling and prediction much easier. 

The fabrication process is electrochemical in nature. When properly biased, 
silicon etches in aqueous HF solution if electronic holes are present in the material. This 
is not the case for n-type silicon in the dark, where no etching takes place. Under hght 
illumination however, electron/hole pairs are generated and etching ^ 

Thus, light intensity is a major parameter to be adjusted to control to the amount of 
silicon etched (see Figures 16-18). If the area ofholegeneratronis 

of the wafer, the pores will grow straight toward it. This is achieved by including the IR 
blocking filter, limiting the depth of light penetration. 

The etching takes place at a certain current density which is a function only of 
temperature and HF solution concentration. This means that if the lattw is kept const *^’ 
the density remains constant Now, by changing the anmunt of hght tomMion 

the etching current does change. This can only mean that the effective ar«t^ ached 
will change since the current density was constant. This area change in turn results 
change in pore diameter. 

But there is a limit to how much the pore diameter can changelf a bare ““P'« 
etched, straight macropores are etched with their location across the w 
distributed Their diameter will center around a size which is determined by the doping 
^SoT and their spring will be demnnined I by the 

oractice one must choose a certain doping concentration for the sample to achieve me 
SpT diameter. Figure 21 shows branching beams, of too much doptng for the 

desired pore size. 

It is also possible to predetermine the pore patent and therdty achieve a much 
more uniform dirtribuuon of pore diameters This pore uutiaion step can ^ done by 
SS? Aching. A phonographic step define, an array ofw.nd ow.m am.sk 
“ched in for example KOH, these result in an of averted pyramds. These 
pyramids arc sufficient to initiate pore growth at those selected sites. 


22 



Thus, to achieve a desired pore pattern, one first chooses an adequate array of 
pore initiation sites. This will determine the pore pitch and array type (linear or 
staggered) of the porous sample. There are limits however in the photolithographic 
process and therefore pore spacing cannot be made arbitrarily small. Then, the actual 
pore diameter is set roughly by the chosen doping concentration of the silicon material 
and more precisely by the amount of light illumination. Since the etching current can be 
easily measured, the pore diameter can then be accurately controlled. Even slight changes 
in hole diameter are possible within the range allowed by the doping constraint. Pore 
depth is matter of time (parabolic function). 


The etching process cannot be made to proceed through the total wafer thickness 
due to two reasons. First, the pores cannot grow towards the source of holes when they 
etch into the generation region; the field confinement of the Ups is no ^onger efifeotve 
Pore branching is the result. Second, the HF etchant fluid would leak out of the setup and 
damage the optics. Instead of etching through the whole wafer the solution employed 
here is to etch the pores through most (~ 90%) of the wafer thickness, and then etch the 
wafer from the backside until the pore region is exposed. 


2. Etching System Setup 

A new etching system was designed and built from scratch. (See Fig. 16 for 
system schematic and Figures 17 and 18 are actual plwtographs >7“™ ;> 

able to use some knowledge and experience gained from our electrochemical «ch*op 
rig, but the inclusion of a light source and other optical parts constituted a new challenge 
since HF as the etching solution does etch glass. 

The rig was designed for operation in the upright position. This allows the 
resuhing^hydrogen* g,, to more eJy detach from the wafer suAce. All ttapm - 
contact 8 with the HF solution are made of Teflon or other compatible mittenal Smce the 
can lead to a temperature nse, a cooling loop was add* (no. shown), 
consisting of tubing, a peristaltic pump and a chiller 

The etching is controlled by a PC and HP Instruments system 
A computer program was written which takes the etching parameters as toe i 
autoZSS/Ss the bias and current output It also include, some satay shutdown 
the system experiences any abnomud opemmg “ndmons. Mly 
note that this system can etch the entire two-inch diameter wafffv S ' S < ^ ^ is in 
eventual implementation of our thermal application into MEMSt *aems _md s 
contrast to the earlier pioneering work done on extremely snudl samples by the Munich 

group [81. 
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Fig. 16: Sketch of etch rig & control loop 



Fig. 18: Photograph of the PC and HP instruments 
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Microfabrication o Porous Structures: Preliminary Results and Success 

a) Preliminary Samples: A couple of lessons had to be learned ot course, before the 
desired pore geometry could be achieved 



Fig. 19: Randomly distributed pores (no initiation) 
showed that the theory practically works. The etch 
depth was fairly predictable and the widening of the 
pores was induced by current-voltage curve 
experiments. 



Fig. 20: Successful pore initiation using a 
mask; without ER blocking filter the pores 
tend to branch. 



Fig. 21: Etched sample showing that too 
high of a doping concentration also leads to 
pore branching of a different kind, 
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b) Success with 5 micron holes and 20 micron pitch array 



Fig. 22: First successful sample etched to a 
depth of 170 microns 



Fig. 23: The above sample after etching the 
remaining silicon from the backside. 



Fig. 24: Even deeper etching reveals some 
pore widening at these etching depths. 
This has not been reported previously and 
can be explained by HF depletion. 


In closing this section we wish to acknowledge Dr. V. Lehmann [9] for his help 
and advice, an in particular, for his suggestions to include a uniform high doping step on 
the backside of the wafer. This small change was a major contribution to the success of 
the porous silicon experiments and resulted in excellent uniformity of the pore geometry 
across the wafer and considerable increase in the illumination efficiency. 


4. Characterization of Porous Wick Properties 

a) porosity and permeability: theoretical calculations 

The porosity of a porous medium, e, is simply defined as the fraction of the total 
volume of the porous matrix that is occupied by the open space [10]. Thus 1-e is the 
fraction of the medium that is occupied by the solid region. The porosity of our given 
porous wick structure shown in Fig. 23 (consisting of 5 micron diameter circular through 
holes arranged in a square in-line pattern spaced 20 microns apart) is thus very easily 
calculated from the geometry and is e = 0.06 Note that the porosity of the porous wick 
structure can be significantly increased by decreasing the spacing between the holes by 
using a different, smaller pitch, hole initiation mask. 

The permeability K of the porous sample is a measure of the flow conductance of 
the fluid-flow through the fluid-saturated porous matrix, and is defined [10-12] according 
to Darcy’s law by the following equation 



dP^ 

dxj 


(15) 


where Ud represents the volume-averaged flow velocity within the porous matrix. Note 
that Ud should not be confused with the actual interstitial pore level velocity u, i.e. the 

local fluid velocity through the open spaces in the porous matrix, but rather Ud is the 
representative porous macroscopic space-averaged velocity through a “sufficiently large” 
volume of the porous matrix containing both the solid and open regions. 

Thus, the permeability, K, can be thought of as the macroscopic manifestation [12] 
of highly viscous flow through the network of pores in the permeable structure, and as 
such its specific value depends on the resulting internal flow arrangement within that 
porous matrix. In general since the flow patterns created within most common types of 
porous medium are highly complex three-dimensional flows, the permeability constant 
appearing in Darcy’s law is usually empirically derived in most situations. 


However, for some very well-defined, simple, and repetitive porous structures of 
known geometry (such as the porous wick structure under investigation as shown in Fig^ 
23) the permeability can also be obtained from the application of the first principles of 
flow of viscous fluids at the pore level (i.e. viscous flow in tubes). Again, it should be 
emphasized that this approach is not feasible for very complex and irregular geometries 
but only can be applied for very simple geometric structures. In short, it can be said that 
no single, overall general relationship exists between the pore structure and permeability, 
the few existing empirical, semiempirical, and first principle-based correlations all have to 
be used within the restrictions for which they have been developed. 
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Fortunately, tor the case on hand, since the porous wick structure is composed of 
a well-defined small-scale network of straight channels, it is therefore possible to 
theoretically derive its permeability as a function of pore size and pore spacing from first 
principle using viscous flow theory To this end the flow through each tube (i.e. pore level 
modeling) is modeled as laminar fully-developed Hagen-Poiseuille flow (The laminar 
fully-developed assumption is reasonable due to the anticipated very small pore Reynold 
number, Rea - 0[1], and due to the rather large L/D = 250/5 ratio.) Recalling from fluid 
mechanics [12], fully-developed laminar flow in a round tube of radius r Q is governed by 


d P _ ' d u 1 d u^j 
dx v<3r 2 r dr J 


(16) 


Solving this equation subject to u = 0 at the cross-section periphery (r = r Q ) yields 


or 


u = 2 U 


\ 2 

r 



ap" 

dxj 


(17a,b) 


where u is the actual local fluid velocity, i.e. u(r), and U is the average or mean fluid 
velocity integrated across the channel (pore) cross-section. 

Now considering once again the entire porous wick structure (global modeling) 
and assuming that Darcy law is applicable, which is a very valid assumption for slow flows 
in porous media, the flow when analyzed as passing through the porous medium as stated 
earlier is governed by the following equation 


K( d?} 

U °~nl dx) 


(18) 


where Ud is the representative porous ('.olume-average) fluid velocity. In terms of the 
average fluid velocity in a single channel. L\ and taking into account the pore density of 
the porous sample, let N/A be the number of pores per unit frontal area (tubes/m^) and 
where A is the total (open and solid) frontal area (m^), the porous fluid velocity can be 
expressed as 
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( 19 ) 


u 


D “ 


N (* r ° ; ) u 

A 


and w hich is simply equal to Ud = eU when one recalls the earlier definition of porosity 
Combining Equations ( 1 7b) - (19) one finds 


N (* r o) 


dP 


A 8^iV dxy 


K 


dP 

dx. 


( 20 ) 


rearranging, permeability is then 


K = 71 r p 4 N = r 0 2 
8 A 8 


( 21 ) 


Equation (21) gives a relation between the permeability and the porosity for our given 
porous wick geometry based on the previous assumptions of fully-developed laminar flow 
and Darcy’s law. 

Figure 25 is a plot of the permeability as a function of the pore size for a given 
pitch pattern based on the above theory As shown, the smaller the pitch, the higher the 
permeability of the porous wick for a given hole size. Smaller pitch means the holes are 
spaced closer together and thus there are a greater number of holes for fluid flow in a 
given cross-sectional area and thus less flow resistance. And, for a given fixed hole pitch 
pattern, we see that the permeability increases with hole size up to a maximum value 
limited by when the individual holes merge, i.e. the pore diameter equals hole pitch for this 
inline arrangement. Here, larger holes again mean more flow area (and higher porosity), 
therefore less flow resistance and thus higher permeability. 

The dotted line in Fig. 25 is the empirical results taken from the latest U S and 
Russian research [13] in sintered power metal wick fabrication and is shown here for the 
sake of comparison. The plotted correlation is a curve fit of actual permeability 
measurements representing the “state-of-the-art” in porous wick fabrication using 
“advanced” fine pore processing technique Although the theoretical permeability of our 
very first successfully micromachined porous sample (solid symbol, 5 micron dia. holes 

with 20 micron pitch; K~ 1.0E-14 m^) is slightly less in magnitude than the best that is 
currently available, this figure clearly re\eals the great potential to easily exceed the 
current limitations by simply micromachining smaller size holes packed closer together 


19 



•«y In**! 



Figure 25 - Theoretical permeability of the micromachined porous wick structure as a 
function of the selected pore radius and pitch with comparison to current 
state-of-the-art wick technology 


b) porosity and permeability: experimental testing 

Porosity of the porous sample are to be determined two ways. Fust, the porosity 
of the porous sample can be easily calculated from its definition given the exact geometry 
of the holes. The representative pore size and pore hole distribution measurements are 
obtained from close up photographs taken under the microscope of the porous sample at 
discrete locations. Circular holes and a repetitive pattern is assumed to extend throughout 
the entire porous specimen and any nonuniformity in the individual pore hole or in the 
pore distribution pattern will result in some inaccuracy in the porosity calculation 
However, judging from the excellent pore uniformity obtained m the successfully 
fabricated porous sample (Fig. 23) we feel that the slight deviations in shape will cause a 
relatively small error in the porosity calculation. Second, it is also possible to the 
determine the porosity of the porous sample by weight analysis by using a very sensitive 
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mass scale if the overall dimensions of the porous sample are known. The advantage of 
this method is that the exact in ternal pore structure does not have to be measured to use 
this technique. However, the uncertainty in the weight measurement and the uncertainty in 
e overall sample s dimensions are the major source of errors to consider when measuring 
porosity using this second methodology. 


The sample B er meability can be determined if the fluid flow through the wick and 
the resulting pressure drop across the wick are known, recall Darcy’s law, Equation (15) 
These quantities are to be experimentally measured using our fluid flow testing loop 
facility shown below in Figure 26. Note that this flow circuit is essentially the same water 
testing loop facility used in the previous microchannel cooling experiments [5] i e the 
same pump, instrumentation, etc., but has been slightly reconfigured with a new test 
section (see insert) specifically designed to securely hold the porous sample in line with the 
water flow. The pressure drop and fluid flow rate are to be measured electronically using a 
pressure transducer and turbine flow meter, respectively, and the data recorded using the 
PC and data acquisition system. We have just successfully completed “leak testing” the 
new test section and are now ready to begin permeability testing of the porous samples 
(See section IV B. 2 for further details.) 



Figure 26 - Porous wick permeability testing loop; insert shows close-up of test section 
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c) effective thermal conductivity of porous silicon wick 

The thermal performance of the wick is paramount to the successful operation of 
the loop heat pipe, and therefore it’s thermal conductivity no doubt will be an important 
factor. In anticipation of this fact we have theoretically estimated, as best we can, the 
thermal conductivity of our porous silicon wick structures by borrowing freely from the 
results of the latest thermal conductivity studies found in the porous media literature [10- 
12, 14-16], 

In general we know [10] that the thermal conductivity of a porous medium 
depends on many variables. Functionally, 


ketr = f[ki, kf, e, microstructure, other considerations) 


( 22 ) 


where 

ketr - effective thermal conductivity of saturated porous medium 
k* - thermal conductivity of solid phase 
kf - thermal conductivity of fluid phase 
e - porosity 


Clearly we are taking a macroscopic point of view (i.e. considering volumetricallv 
averaged or global quantity, not pore level or local value of thermal conductivity) in 
dealing with keff . (In treating the wick structure as a porous media the solid phase is 
silicon and the open pore region is the fluid region, liquid or vapor.) 

There are two extremely simple 1-D models which establish the upper and lower 
bounds for keff of a porous medium (in general) and are based upon assuming either a pure 
parallel type connections (Eq. 23) or a pure series type connections (Eq. 24) of the 
constituents making up the porous medium. Basically this amounts to assuming the porous 
material is a 1-D heterogeneous material of known structure. 


Temperature Ti 

Tqnpggggt T| 

k S 



k. 

kf 


kr 


Temperature Tj Temperature Tj 


(Parallel) 


(Series) 


keff =(1 -e)k, +ek f (23) 


1 (1-S) t 6 

keff k, + k f 


(24) 
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Often, however, porous matenaJ is composed of packed beds of particles and as 
sue contains a combination of both types of heat transfer paths (series and parallel) 
Empirically this arrangement has been shown to behave according to P 


L- _ L.O-e) \. £ 

K eff “ K f 


(geometric mean) 


(25) 


as long as kr is not two much different than k, , with the value of kc***™*™ falling 

m between the ^ lumts J ust described. Also note that regularly or randomly 
packed beds to a very good approximation can be treated as a homogeneous medium with 
respect to the effective thermal conductivity, i.e. no directional dependence. 


More accurate formulas for kdr are available in the literature but become 
increasingly more complex, requires additional information, and are more “system 
specific” applicable to a particular type of porous medium. For example better 
correlations exists (based on 2-D and 3-D models) specific for a given particle geometry 
(spheres, rods, etc.) with known packing arrangements (simple cubic, body center cubic 
etc.). Also, if packing is random but tend to take on a certain distribution, let’s say four 
different lands of packing structures, a rang? of lwrcan be theoretically predicted based on 
the probability distribution of finding those certain types of packing structures within the 
porous medium. Equations (26)-(27) listed below is of this last type and gives the upper 
and lower bounds for the traverse conductivity kest™*. for a two-dimensional anisotropic 
fibrous material composed of needle-shaped fibers, all being arranged in parallel with each 
other, in a gaseous or solid matrix of different thermal conductivity [14]. 


VMr 


1+MP-O 

, (5,(1-6 V XP-1) 2 ) ]' 


[ ^ 1+ 5 v (0-l)ll+(|i-l)6 v +2O-lXl-25 v )G]j 
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2(1 - p)8 v + 2P + 4< p - IX28 V - 1)G 


-i-i 


(Upper limit) (26) 


(Lower limit) (27) 


Where 

0 = ~ G=0.25 for Circle shaped cells 

* r 

Note that the recommended effective thermal conductivity (according to reference [14]) in 
the axial direction (parallel to the fibers) for the same anisotropic fibrous porous medium 
structure is simply given by the earlier Eq (23) Hence, to summarize, for this stated 
fibrous arrangement it’s concluded [14] that the very simple parallel 1-D model is 
reasonable for predicting kegMui but the kdfnv« values range somewhere between the 
upper and lower 2-D limits as given by Eqs ( 26 ) and (27), respectively (which are 
significantly larger in magnitude than the I D senes model estimate of Eq. <24) due to 
multi-dimensional heat conduction in the medium ) 



These different thermal conductivity models just described are used to estimate the 
effective thermal conductivity of our porous silicon wick structure, saturated with water in 
Table 1, and saturated with water-vapor in Table 2, for select values of porosity The 
results of the wick saturated with water is also shown plotted in Figure 27 Inspection of 
Tables 1-2 along with Figure 27 reveals the following: 

1 > rhe value of keff of the wick structure saturated with water or saturated with water- 
vapor is approximately the same for k^.** (Eq. 23), keg*,™. up p« lunit (Eq. 26), and 
keff, geometnc mean (Eq. 25) models, i.e. compare the values of Table 1 and Table 2, but 
differ significantly for keg travene low limn (Eq. 27) and kegam, (Eq. 24) models. 

2) keg exiai (Eq. 23) varies linearly with porosity (thus is a strong function of e). 

3) keg travene 2-D upper and lower limits, Eqs.(26) - (27) » k«gi-D am, Eq. (24). 
(Hence the 1-D series model greatly underestimate the value of kagtraww). Also note 
that kegtraveraa dependence on porosity diminishes as e increases, now dominated by the 
value of the low conductivity fluid). 

4) keggeomatric mean, Eq. (25), is always within the upper and lower bounds of kagtrmvena 
and is a much simpler equation to use to calculate the effective thermal conductivity. 

5) kegaxial > keg traveree 2-D upper limit for all Values of porosity. 

In applying the above results to the wick geometry there are two additional facts that one 
should keep in mind: 

6) The wick geometry is not exactly the same as that of the 2-D fibrous insulation 
system modeled and analyzed in ref. [14]; however both systems are reasonably similar 
in shape to one another to a certain degree, note also that ref. [14] geometry is the 
closest that can be found to the wick structure in the literature without doing a more 
extensive search. 

7) The results of reference [14] is based solely on a theoretical study only and was not 
validated experimentally at the time (although latter reference to this study imply that 
the reported results are reasonable.) 

Hence, given these last two comments imply that one should use caution in applying the 
results of Eq. (23) and Eqs.(26)-(27) (or Eqs (23) and (25)) to determine and 

kegtnwera*, respectively, of the porous silicon wick. The results obtained using the above 
equations therefore should be viewed only as rough estimates only, not actual values of 
the effective thermal conductivity of the porous wick structure. Finite element modeling 
[16] are recommended (and/or actual thermal conductivity testing of the wick is 
suggested) if more precise values of thermal conductivities of the wick structure are 
needed. (However we do feel that the equations listed above reveal the proper trends, i.e. 
effect of porosity, thermal conductivity of the saturating fluid (water vs. water- vapor) 
directional dependence (axial vs. traverse), etc , and, in fact, can provide reasonable first 
cut estimates.) 
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Table 1 

K.-128 W/m *C 

Kr*0.648 W/m 

°c 

(Silicon / Water) 

Porosity 

Eq. 23 

CM 

cr 

UJ 

Eq.25 

Eq. 26 

Eq. 27 

0.2 

102.530 

3.176 

44.472 

82.359 

13.946 

C.5 

84 314 

1.289 

9.107 

32.807 

2.528 

* 7 

W . r 

38.840 

0.924 

3.164 

12.379 

1.309 

0.9 

13.360 

0.720 

1.099 

2.037 

0.798 

Table 2 

K.»12tWAn*C 

Kv24.8x10 J W/m*C 

(Silicon / Water 

■Vapof) 

Porosity 

Eq. 23 

Eq. 24 

Eq. 25 

Eq. 26 

Eq. 27 

0.2 

102.405 

0.123 

23.118 

81.937 

0.811 

0.5 

64.012 

0.049 

1.774 

32.031 

0.098 

0.7 

38.417 

0.035 

0.320 

11.553 

0.050 

0.9 

12.822 

0.027 

0.058 

1309 

0.030 



Fig. 27 Thermal conductivity (theoretical) of the micromachined silicon porous wick 
structure saturated with water as a function of porosity 





IV. PLANNED W ORK fYear 3 Goals! 


B. Porous Wick Structure 

1. Microfabrication of Porous Silicon W ick Structures 

Our goal next year is to fabricate a second set of porous silicon test specimens, 
again with straight through holes, but this time the pores fabricated in the silicon wafer 
will be smaller and spaced much closer together In particular, additional pore initiation 
masks will be fabricated which allow for pore arrays with a 4-micron pitch. This is at the 
lower limit of our photolithography equipment Both linear and staggered pore arrays will 
be made and the pore diameter will be reduced to about 1 micron in these samples. 

The slight widening of the pores as observed in our earlier results (Fig. 24) will 
also be investigated and studied We will control this behavior in our earlier 5-micron 
diameter pores samples by adjusting the light intensity. Then, if necessary, it will also be 
dealt with accordingly in the latter samples containing the smaller hole arrangements. 
Since this effect has not been previously reported in the literature our findings may result 
in a short publication. 

The fabrication of a porous silicon wafer with the much finer pore patterns will 
essentially conclude the macro porous silicon fabrication efforts here at UC, although it 
needs to be stated that much more work still remains to be done before a micromachined 
silicon porous wick can be successfully integrated into a loop heat pipe. (See section B.3 
for further discussion.) 

2. Testing of Porous Silicon Samples 

The porosity and permeability of the porous silicon samples with be measured here 
at UC using the techniques and test set-up as described in section III.B. 4.b. Testing of the 
5 micron pore sample is currently in progress; the smaller pore samples will also be 
similarly tested when available. It is our hope that the measured permeability will closely 
follow the theory but any major deviation in the data from theory will be resolved 
Thermacore and Texas A&M are expected to confirm our results by doing repetitive 
testing on similar samples (fabricated by UC) using their own testing facilities. The 
effective thermal conductivity of the samples will not be measured. 

3. Microfabrication of Advanced Porous Structures 

As mentioned in the Project Description section, there may be some advantage to 
varying the pore diameter across the thickness of the wick. This should allow one to 
operate under normal conditions with one radius of curvature, and utilize another radius of 
curvature to prevent the meniscus from receding into the liquid. Also, by shortening the 
length of the smaller diameter section, the hydrostatic pressure drop can be reduced 

Several discussions have occurred between the University of Cincinnati, NASA 
Lewis, DTX and Texas A&M Two or three concepts have been proposed. Before any 
final design, however, DTX and Texas A&M will have to perform some detailed 
calculations to specify the actual dimensions which are desired. The University of 
Cincinnati can then determine -if it is possible to fabricate. Of course several iterations may 
be necessary before a workable design can be reached 
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V. CLOSING 


Although many multi-chip module (MCM) approaches have been tned or used in 
recent years, the one herein proposed is quite revolutionary, but yet realistic Work is 
progressing at a steady rate on schedule and major achievements have already been 
realized at this point in time Ultimately, the micromachined integrated MCM thermal 
management system should demonstrate that the chip temperatures can be maintained at 
the present levels while increasing the volume density by as much as two orders of 
magnitude at the chip level. The resulting savings in mass and volume we feel will prove 
to be extremely beneficial to NASA. 

Also, an additional endeavor that is to be accomplished during the course of this 
project is to use MEMS technology to fabricate "porous" silicon, with specific 
"micromachined" geometric characteristics, to form the ultimate in wick structures for use 
in heat loops. To the best of our knowledge porous silicon has never before been 
considered for such a thermal application but has rather been investigated for its 
luminescent properties (not fully understood). However, since the porosity, anisotropy 
and pattern configuration can be controlled it could revolutionize the area of 
semiconductor cooling, thermosyphon evaporation, etc. Moreover, the planar technology 
is a batch-processing technique, which will ultimately be very inexpensive. 

Of course successful creation of porous silicon wicks will be of great benefit for 
utilizing this MCM concept in space applications. Eliminating the pump will further 
reduce the mass and volume of the complete thermal management system. 

m 

Thus, we feel that this project embodies several rather revolutionary and potentially 
very pragmatic areas of thermal development by cross-fertilization of several disciplinary 
areas heretofore overlooked. 
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Agenda 


1. MCM - Thermal Management System 

* Microfabrication 

board etching, chip separation, sensors, etc. 

* Thermal Analysis 

contact resistance, 3-D thermal model 

\ 

2. Porous Wick Structure 

* Microfabrication 

* Permeability Testing 



University of Cincinnati 




Porous Wick Structure 


□ microfabrication facility 

u 0 5pm x 20pm pitch fabrication 

□ new results 

□ permeability testing 
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Porous Wick Structure 
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Porous Wick Structure 
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Fabrication Sequence 


1 -2. Oxidation. BOE on Backside 7. Pore Initiation Etching (KOH) 

Si0 2 
n - Si 


3. N* - Diffusion 


8. Porous Silicon Etching (HF) 
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4 Metal Deposition 


9., 10. Al Etching, Reoxidation 
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n* - Si 
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5. Metal Lithography and Etching 


11., 12. Backside BOE, Wafer Thinning (TMAH) 




6. Pore Initiation Lithography 


13. BOE 
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Porous Wick Structure 





University of Cincinnati 



First successful sample etched to a depth of 170 microns 


Deeper etching (260 pm) reveals some pore widening. 



The above sample after etching the remaining Deep etching with adjusted etching current yields 

silicon from the backside straight pores 
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Porous Wick Structure 
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New Results 


□ 0 1-3|jm x 4pm pitch fabrication 

□ new mask, parameters 



0 2.5pm x 4pm pitch 0 1.1pm x 4pm pitch 
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Porous Wick Structure 
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Permeability Testing 



Modified fluidic testing setup for 
permeability measurements. 

Porous sample test section, assembled 
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Porous Wick Structure 
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Permeability Testing 
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Porous Wick Structure 






Theoretical Permeability Development 


Pore level modeling 

• Fully-Developed Laminar Flow in a Round Tube ( Flagen-Poiseuillie flow): 


Governing Equation 


dP_ 

d x 


S 2 u 1 du 
\dr~ r dr) 


( 1 ) 


Solution 


U = 


r 2 ( dP\ 


8 ju v dxJ 


( 2 ) 


Global modeling 

• Darcy model flow in the porous wick 

Kf dP\ 
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( 3 ) 
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■ v (^r) 


U = eU 


Where 


(4) 


N/A the number of pores per unit frontal area (tubes/m 2 ) 

A the total frontal area (m 2 ) 

s the porosity 


Combining (2) & (4) 
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Permeability Testing and Data Reduction 
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Permeability, Kim 2 ] 


Permeability Testing 
(Experimental Versus Theory) 
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Permeability Variation with Pore Radius 



Summary 


• rr.icrochannels and milli-slot have been 
fabricated, tested and thermal model verified 

• prototype MCM module with micromachined 
pockets and chips has been fabricated 

• 3-D thermal model of MCM module is under 
development 

• contact resistance experiments, 3-D chips-to- 
substrate, are almost ready to begin 

• 5 pm holes x 20 pm (pitch) silicon porous 
wicks have been successfully fabricated 

• permeability measurements of porous wick 
test specimens are ongoing 

• 1 pm holes x 4 pm (pitch) silicon porous 
wicks are under fabrication 


• future work 
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ABSTRACT 

This paper reports a new concept of building a coherent porous wick structure in planar silicon wafers for capillary pumping 
in Loop Heat Pipes (LHPs). By utilizing Macro Porous Silicon (MPS) fabrication technology, with pore sizes in the order of 
a micron, precise control of the pore dimensions of the wick is possible. In addition, by using MEMS (Micro Electro 
Mechanical Systems) fabrication technology, LHPs can ultimately be integrated into electronic packaging systems, or indeed 
into the silicon microchips themselves. The MPS samples were fabricated and tested for capillary pressure and permeability. 
The test results closely matched the analytical models that were derived from basic physics. 

Keywords: Macro Porous Silicon (MPS), Loop Heat Pipe (LHP), wick, electronic cooling, microchannel flow 

1. INTRODUCTION 

Enhanced cooling of present and next generation electronic packaging systems is crucial. The realization of increased speed 
and higher device density result in ever-higher heat fluxes, calling for design of fully integrated packaging-cooling systems. 
This work reports a portion of a larger ongoing effort to develop a unique all-silicon MCM (Multi Chip Module) design [ 1 ] 
which includes several concurrent cooling schemes, such as forced liquid cooling in ultra deep microchannels [2], passive 
cooling with Loop Heat Pipes (LHPs) and more traditional air fan cooling. The passive nature of the heat pipe operation 
renders them particularly attractive, since no external driver is required. The porous wick in the device is itself the pump of 
the LHP, using only the very heat it removes. 

First, this paper covers heat pipe and LHP fundamentals and the particular importance of the wick. Then, the characteristic 
parameters of capillary pressure and permeability are introduced and applied to the coherent wick. Macropore fabrication is 
included with particular attention paid to the necessary and important fabrication details. The testing results are then 
compared to the theoretical predictions. Finally, an integrated evaporator concept for MCMs is introduced. 

2. COHERENT WICK DESIGN 


Motivation 


A heat pipe is a sealed device in which a working 
fluid circulates and constantly changes from the 
liquid to the gaseous phase, and vice versa 
(Figure 1). The evaporated liquid flows from the 
evaporator to the condenser end by the pressure 
drop created from the evaporation and 

condensation process. The liquid is returned to 
the evaporator by the capillary pressure in the 
wick. The input heat energy is thus converted 
into latent heat in the evaporator, moved in the 
gaseous phase to condenser, and then converted 
back into heat energy by the evaporation process, 
and the heat is then ejected at the condenser end. 
The evaporation and condensation process takes 
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Figure 1 . The basic cycle of a heat pipe. 



place at almost the same temperature, i.e. it is virtually isothermal. In a 2D analogy, one often compares a heat pipe with a 
highly thermally conductive material. The finite heat flux q is transported at a small temperature difference AT . resulting in a 
high effective thermal conductance G fhM y: 
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Heat pipes offer a passive way of rapidly transporting energy away from a heat source. The term “passive" is used here in the 
sense that no external action, such as pumping, is required for operation. The heat itself pumps these devices. This has certain 
advantages over other heat removal techniques. No extra devices (such as pumps or other actuators) are needed, and no 
additional power needs to be introduced as heat into the system. 


The distinction between an ordinary heat pipe (HP), where a wick continuously lines the inner surface of a single pipe, and a 
capillary pumped loop (CPL) is that in a CPL interconnecting plumbing separates the condenser and evaporator (including 
the wick). The length of an ordinary heat pipe is limited, because the whole liquid return path is through the wick, which 
results in high viscous losses. The capillary pumped loop overcomes this shortcoming, but it is not inherently self- starting. 


A recently developed type of heat pipe is the Loop Heat Pipe (LHP) [3]. It is similar 
in design to a Capillary Pumped Loop because it uses separate tubing for the liquid 
and vapor flow paths, which allows the thermal energy to be moved over relatively 
large distances (Figure 2). The wick is located at the evaporator and it serves as the 
engine, which achieves continuous circulation. The LHP is unique because the 
overall design provides for a reliable start-up. Working LHPs have been made using 
cylindrical and planar evaporator configurations. 

The porous wick is without question the most important part of a LHP. It transfers 
the energy to the liquid-vapor interface ("meniscus”) and maintains separation 
between the vapor and the liquid by capillary forces. The porous wick thus 
determines the pumping capability of a LHP like no other part. MEMS (Micro 
Electro Mechanical Systems) technology opens up new possibilities in wick design, 
as well as modeling. Traditional, sponge-like porous wicks have a randomly 
distributed pore size, and the flow path is irregular. Consequently, the mathematical 
description is statistical and semi-empirical. On the other hand, the coherent nature 
of wicks fabricated by the Macro Porous Silicon (MPS) etching technique allows a 
relatively simple analytical description of the parameters. 
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Figure 2. Schematic of a Loop Heat 
Pipe, redrawn from [3] 


The capillary action pumps the liquid from the condenser to the evaporator. The capillary pressure has to overcome all other 
pressure losses in the system. Wetting and surface tension result in a curved liquid surface, or meniscus. Macroscopically, the 
pressure drop A P t across the meniscus is described by the surface tension <r, the wetting angle 0, and the pore radius r p [6]: 


A P = 


2 o cos0 


( 2 ) 


Practically, one often introduces an effective surface tension, /rcos0, since it is difficult to separate the two components in 
boiling experiments. In order to achieve a high capillary pressure, a small pore radius, good wetting and a liquid with large 
surface tension are required. 


Permeability 

The permeability K is a measure of the liquid flow conductance through the porous wick. It is used in the expression 



( 3 ) 



where u n is the macroscopic, volume-averaged flow velocity. // the dynamic viscosity of the liquid, and cP/cx is the axial 
pressure gradient. Assuming fully-developed laminar How in the pores ot the wick one can show [4] that, for an array ot 

straight pores in a wick of porosity e (the ratio of pore volume to total wick volume), the permeability is. 

■> 

K = ^£. < 4 > 

8 

Thus, to maximize the liquid flow through the wick for a given pressure drop A P, one must choose a liquid of low viscosity, 
minimize the wick thickness d wick =Ax, and maximize the pore radius and porosity. Note that the latter requirement contradicts 
that for a high capillary pressure in Eq. (2). 


Fabrication 


MPS technology was first developed by Lehmann [7]. 

It allows the formation of coherent arrays of straight 
pores in n-type silicon, with a pore diameter ranging 
from 0.2-20|im. 

The fabrication process is electrochemical in nature. 

When properly biased, silicon etches in aqueous HF 
solution if electronic holes are present in the material. 

This is not the case for n-type silicon in the dark, thus 
practically no etching takes place. Under illumination 
electron/hole pairs are generated and n-type silicon 
does etch. Therefore, light is one major control 
mechanism for the amount of silicon etched (see 
Figure 3). The primary directions of pore growth are of 
type <100>. If the area of hole generation is limited to 
the backside of the wafer, the pores tend to grow 
toward the source of holes. The pore tips contine the 
electric field and consume all the holes for etching. 

Such a constraint of the generation region is achieved 
by including the IR blocking filter, limiting the light 
penetration depth to a small thickness on the substrate 
backside. On a (100) wafer, the vertical pore growth 
direction is one of the <100> directions. Therefore, 
deep and coherent MPS can only be formed, at this point, on n-type (100) silicon wafers. Other researchers were able to form 
MPS on p-type wafers without illumination [8,9], but the resulting pores are not as deep and uniform as the ones described 
here. 

The etching at the pore tips takes place at a certain current density, J PS * which is only a function of temperature and HF 
solution concentration (see Appendix). This means that if the latter are kept constant, the current density remains constant. 
Now. by changing the illumination intensity, the etching current does change. This must result in a change of the effective 
area that is etched, since the current density is constant. This area change translates into a change in pore diameter. Yet there 
is a limit of how much the pore diameter can change. The range of allowed pore sizes is determined by the substrate doping. 

There is also a way to predetermine the pore pattern and thereby to achieve a very uniform pore array. Photolithography and 
subsequent anisotropic etching can be used to achieve a pore initiation step. The resulting pyramidal etch pits [10] are 
sufficient to concentrate the field at their tips and initiate pore growth at these pre-selected sites. 

In summary, to achieve a desired pore pattern, one first chooses an adequate array of etch pits for pore initiation sites This 
will determine the pore pitch. Then, the actual pore diameter is set roughly by the chosen doping concentration ot the silicon 
material and '"fine-tuned" by the amount of illumination. Since the etching current can be easily measured, the diameter can 



Figure 3. Schematic of the macropore etching setup 




he accurately controlled. Small changes in pore diameter are possible within the range allowed by the doping constraint. 
Finally, the pore tip etches with an etch rate that is proportional to the current density J r s- Thus, the desired pore depth is 
achieved by timing the etching process. 

The etching, however, cannot quite continue through the whole wafer thickness because ot two reasons: (1) The pores at 
some point cannot grow towards the source of holes as they etch into the generation region, because the field confinement of 
the pore tips is not effective any longer. Thus, as the pores approach the broad backside field plane, they tend to branch. (2) 
The HF would leak out of the setup and damage the optics in our system. Instead of etching pores through the entire water 
thickness, the pores are etched through most of the wafer, and then, the wafer is selectively etched from the backside using 
TMAFI or another appropriate etchant, until the pore region is reached. 


Lehmann [7] already pointed out that the etching process leads to an effective decrease in HF solution concentration at the 
pore tips, resulting in a parabolic development of the pore depth over time. It was found in the present study that this also 
leads to an increase in pore diameter with increasing depth if the etching current is kept constant (Figure 4). Based on the 
parabolic development of the pore depth, the concentration change can be calculated and the current can be appropriately 
adjusted (see Appendix). Samples etched with these improved parameters exhibit a constant pore cross-section (Figure 5). 



Figure 4. Cross-sectional view of a sample etched at a 
constant current. The diameter increases with 
depth. 



Figure 5. A sample etched with a current that was 

adjusted to yield a constant pore cross-section 
(appendix). 


Thus far, wicks w-ith pores of about 5-7pm pore diameter and a pitch (center-to-center distance) of 20pm were fabricated. At 
the moment, our effort is centered towards fabrication of pores with a smaller diameter ( 1 pm) and a smaller pitch (4pm). 




Testing setup for Measurement of Permeability and Capillary Pressure 

Deionized water was used for the flow 
tests, where the pressure drop across the porous 
sample was measured using a manometer, and 
the How rate was determined by collecting the 
volume of liquid over a certain time duration 
(Figure 6). 

The capillary pressure was obtained using the 
same setup. Initially, with only liquid on the 
pump side and air on the other side of the porous 
wick sample, menisci formed at the interface. 

The liquid initially did not flow, and pressure 
builds up. Only when the pressure exceeds the 
capillary pressure does the liquid begin to flow, 
and the pressure drop across the sample 
decreases suddenly. The pressure falls to the 
value determined by the permeability of the 
sample. The maximum pressure obtained before 
the flow starts is the capillary pressure. 

3. RESULTS 

As Eq. (3) suggests, the permeability of the sample can be deducted from the slope in the flow / pressure drop diagram. For 
one test run, a line was fitted into the data field (Figure 7) and the permeability calculated. The results of the tests along with 
the theoretical predictions, as well as test data from conventional sintered porous metal wicks [4], are plotted in Figure 8. 
Although permeability data were collected from samples of only one wafer, the results so far confirm the theoretical 
predictions. Even with its low porosity (7%), in terms of permeability the coherent MPS wick performs as well as the 
conventional porous metal wicks (themselves a fairly recent development) which frequently have a porosity value in excess 

of 50%. 

The capillary pressure tests conducted to date are not very precise since the tests were primarily conducted to collect flow 
data for permeability measurements (TABLE 1). The additional data shown in the table are measurements taken by 
Thermacore [11] (commercial developer of porous metal wicks) with the so-called bubble-point test. Their results also show 
a rather high degree of uncertainty. However, the theoretically predicted radii derived from the measured capillary pressure 
are also close the SEM measurements. 



Figure 6. Testing setup for capillary pressure and permeability. 
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Figure 7. Example of a flow rate versus pressure drop test 
The slope is inversely proportional to the 
permeability. 
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Injure x Permeability vs. pore radius for two different pore 
initiation patterns (masks), as well as conventional 
porous metal wicks. The theoretical and test results 
match well. 




TABLE 1. 

Capillary pressure test results 


Test Lab 

Wafer 

(qm) 

(from SEM) 

fluid 

. 

Test# A/Hpsi) 

r.„,„(pm) 
(from A p) 

Univ. of 
Cincinnati 

PS20 

3.25 

water 


6-8 

2.6- 3.5 

l — 

Thermae ore 

[Ml 

PS 17 

2.5 

acetone 

5 

1.17 

2.825 ‘ 

6 

1.27 

2.6 

7 

1.55 

2.225 

8 

1.41 

2.44 


4. INTEGRATION 


The coherent MPS wick must be ultimately integrated into a working LHP. This would be necessary to clearly demonstrate 
the feasibility of using a MPS wick. 


Anisotropic etching and wafer bonding techniques can be used to create a flat evaporator (Figure 9). On the vapor side 
energy is conducted to the wick by an array of mesas from the hot plate. This plate is anisotropically etched from a ( 
silicon wafer using KOH and appropriate comer compensation techniques. The volume between the mesas provides the 
vapor flow path. This hot plate could be attached directly to an MCM, or be itself part of the all-silicon MCM [ 1 ]. The porous 
silicon plate and the bottom plate form part of the liquid reservoir. 
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Figure 9. Cross-sectional view of an evaporator fabricated by MEMS 
technology. 


This integrated planar evaporator will replace the conventional, cylindrical 
evaporator of a working LHP. It is extremely important that the evaporator is 
flat for MCM cooling applications, and so far conventional porous metal 
wicks could not achieve this goal because only cylindrical evaporator sections 
have been made. Ultimately, in the future it is envisioned that a whole LHP is 
fabricated by MEMS (Micro Electro Mechanical Systems) technology. 


5. CONCLUSIONS 

MPS was used to create a new type of wick for use in LHPs. The good 
dimensional control ot pore size and pattern offered via state-of-the-art 
MEMS-based micromachining allows optimal fine-tuning of the wick’s 
parameters, such as permeability and capillary pressure. The feasibility of 
fabricating a novel wick by MPS etching was demonstrated. Relative simple 
analytical modeling of important performance parameters of the wick was 
carried out and was verified experimentally. 




Operation of the mieromachined silicon wick in an actual LHP has yet to be confirmed. A parametric study is also suggested 
to optimize the wick performance, as well as the thermal and mechanical integration with a LHP. 

In terms of permeability and capillary pressure, the porous silicon wick is at least as good as are conventional porous metal 
wicks, even at a far lower porosity. Being a MEMS-type component, the porous silicon wick has several advantages. It is 
inherently planar, can be made very thin, and it can be integrated directly and intimately with other MEMS components. 
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8. APPENDIX 


A: coefficient 

c: HF solution concentration in wt.% 

c 0 : initial HF solution concentration 

C: pre-exponential factor (3300 A/cm 2 wt.% 273 ) 

d\ pore cross-sectional diameter 

e\ elementary charge 

E a \ activation energy (0.345eV ) 

/: etching current 

J rs \ critical current density 

J rso \ initial critical current density 

k: Boltzmann constant 

L : pore length (depth) 

n: dissolution valence (2.4) 

N Sl : atomic density of silicon 



r. ume 

T, tempera lure in K 

r: pore growth velocity (etch rate) 

r„: initial pore growth velocity 

Lehmann |7) determined the Arrhenius-type critical current density (at which the silicon is etched at the pore tips): 


J rs = CV 2 exp 



and the etching velocity, or etch rate 


( 5 ) 


n-(-e)N Sl 

Lehmann also experimentally determined the following description of the pore depth with time, for a 2.5 wt.% HF solution at 
room temperature: 

t _ 1.89 ,, 0.0029 , 2 (7) 

— — l_j ~ ^ Lj . 

min pm pm" 


Using Eqs. (5) and (6), one obtains the following relation for the decrease in HF concentration at the pore tips: 



This reduction is significantly different from the linear fit of Lehmann, for an etching time more than about 6 hours or 1 50pm 
deep (Figure II)! Now, Eq. (5) yields the relation for correcting the current: 



Also the pore widening such as in Figure 4 (when etching without this current correction) can be predicted (Figure 1 2): 



For different HF concentrations and temperatures a relationship similar to Eq. (7) has to be found. The general form is 

t = v n ~'L + A(c,T)L 2 

Thus, a single test is sufficient to determine the coefficient A. and it can be decided if a current correction is necessary. 
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Figure 1 1 Decrease in HF concentration after Lehmann 
and with the present newly refined model. 



Figure 12. Calculated relative increase in pore diameter 
when etching with a constant current. 



